Anatomy of hypothalamic and diencephalic nuclei involved in seasonal fertility regulation in ewes

In this study, we describe in detail the anatomy of nuclei involved in seasonal fertility regulation (SFR) in ewes. For this purpose, the intergeniculate leaflet of the visual thalamus, the caudal hypothalamic arcuate nucleus, and suprachiasmatic, paraventricular and supraoptic nuclei of the rostral hypothalamus were morphometrically and qualitatively analyzed in Nissl-stained serial sections, in the three anatomical planes. In addition, data were collected on calcium-binding proteins and cell phenotypes after immunostaining alternate serial sections for calretinin, parvalbumin and calbindin. For a complete neuroanatomical study, glial architecture was assessed by immunostaining and analyzing alternate sections for glial fibrillary acidic protein (GFAP) and ionized calcium-binding adapter molecule 1 (IBA1). The results showed a strong microglial and astroglia reaction around the hypothalamic nuclei of interest and around the whole 3rd ventricle of the ewe brain. Moreover, we correlated cytoarchitectonic coordinates of panoramic serial sections with their macroscopic localization and extension in midline sagittal-sectioned whole brain to provide guidelines for microdissecting nuclei involved in SFR.


. Introduction
Sheep are ruminants and seasonal breeders. As a result, sheep farming is misaligned with market demands (e.g., continuous production throughout the year), and the income derived from these seasonal animals is lower than that derived from non-seasonal animals (1). For this reason, one of the most important goals for seasonal animal husbandry (e.g., sheep) is to achieve continuous production throughout the year by shortening the periods between lambing, thereby improving production rates (2, 3).
Seasonality, an adaptive process whereby animals adapt their physiological functions (e.g., reproduction) to the environmental conditions (4), has been studied in several species ranging from invertebrates to vertebrates (5)(6)(7). In particular, research has focused on understanding the neurological basis of seasonal breeding among farm animals (e.g., sheep and goats) (8)(9)(10). Thus, structural plasticity changes, cell proliferation and migration in the sheep hypothalamus indicate that seasonal light changes induce a deep neural reorganization .
/fvets. . (11)(12)(13). Furthermore, recent studies have shown that polymorphisms in genes encoding melatonin 1A and 1B receptors (MTNR1A and MTNR1B, respectively) are involved in reproductive seasonality (14)(15)(16)(17). Concurrently, anatomical and physiological studies have implicated melatonin in the control of seasonal breeding or estrous cyclicity (18,19). Based on these findings, attempts have been made to regulate the circadian cycle of sheep by changing their natural photoperiod in combination with melatonin implants and hormone treatments (20)(21)(22). In immunohistochemical studies, kisspeptin, neurokinin A and dynorphin B, called KNDy neurons, and Gonadotropin Releasing hormone (GnRH) neurons have also been identified in the sheep hypothalamus, highlighting their key role in the endocrine regulation of seasonal breeding (8,(23)(24)(25). However, analyzing interactions and the complex neuronal machinery involved in the hypothalamic regulation of seasonality requires understanding its neuroanatomical underpinnings. The main nuclei involved in circadian rhythm regulation in rodents and non-human primates are the hypothalamic suprachiasmatic nucleus (SCh) and the diencephalic intergeniculate leaflet (IGL) nucleus of the thalamus (26). Other hypothalamic nuclei also involved in seasonality include the supraoptic (SO), paraventricular (PVN) and arcuate (ARC) nuclei (27, 28). With slight differences in morphology and spacial distribution, these nuclei are present in all higher mammals, including sheep (11,26,(29)(30)(31). Full neuroanatomical and neurochemical studies of hypothalamic nuclei involved in seasonal regulation (e.g., circadian rhythms and seasonal changes in GnRH regulation) must be conducted to gain further insights into the neural basis of fertility regulation in sheep (24, 32, 33).
A few anatomical atlases of the ruminant brain (e.g., sheep, goat, and alpaca) have been previously published (34)(35)(36)(37). These studies have used a wide range of methods, including gross anatomy, high-resolution 3D Magnetic Resonance Imaging, and immunohistochemistry or Nissl-stained sections. A highly detailed general description of sheep diencephalic nuclei was reported in a doctoral dissertation dating back to 1976, using coronal thioninestained seriated sections (38) however it was not interpreted in the context of the neuromeric model. Thus, further research should be conducted to analyze in detail the sheep neuroanatomical structures involved in seasonal breeding.
In this article, we describe the anatomical analysis of SFR nuclei using Nissl-stained serial sections and immunocytochemistry for CBPs. For this purpose, we selected three CBPs, namely Cr, Cb and PV, which have been widely used as markers for cytoarchitectural analysis and for neuronal phenotyping in different vertebrate species (39,40).
Considering the above, the present study aimed at, (1) determining the morphological features and relative position (coordinates) of seasonal fertility regulation (SFR) nuclei in the ewe diencephalon and hypothalamus by analyzing Nissl-stained serial sections in the three anatomical planes based on high-resolution panoramic photographs, (2) studying the immunohistochemical calbindin (Cb), calretinin (Cr) and parvalbumin (Pv) profile of SFR nuclei, (3) analyzing the morphology, distribution, and cell density gradients (glial architecture) of microglial and astroglia cells in the ewe hypothalamus, and (4) generating a macroscopic map with the specific coordinates and extension of SFR nuclei to enable their accurate localization for microdissections.
. Materials and methods

. . Experimental groups
Six adults female Churra sheep (age = 16.60 ± 0.25 months); mean live weight [(LW) = 70 ± 6 kg] housed under natural light, temperature, and humidity conditions, fed ad libitum and provided with unrestricted access to clean drinking water, were used in this study over farm natural conditions. The experiments were performed in accordance with Spanish (RD 53/2013) and European (63/2010/EU) directives, and animal maintenance and care met all requirements of the agreement on the use of animals in scientific experimentation in Spain (Confederation of Scientific Societies of Spain -COSCE). Accordingly, the animals were sacrificed at the slaughterhouse MACRISA, in Medina de Rioseco, Valladolid, Spain, in April 2021, by neck cutting, severing the carotid and jugular veins, without prior drug administration. Heparinized 5-ml tubes were used to collect blood samples by jugular venipuncture; samples were immediately centrifugated at 3000 × g for 20 min, and plasma was stored at −20 • C. Hormonal levels of progesterone was assayed by radioimmunoassay (RIA).All samples were run in a single assay. In all specimens, the progesterone levels were homogeneous, ranging from 15 to 20 mg/l.

. . Animal samples
After removing the skin, a horizontal cut was made in the skull from the external occipital protuberance to the upper edges of the orbits using an ultra-fast surgical saw with a 5 cm diameter radial disc. Using this approach, the dorsal superficial area of the cortex was cleanly cut, thus showing in depth the ventricles in the middle line maintaining intact the diencephalic area. This procedure provided access to the ventricles for a faster and more effective fixation. After tilting the skull backwards with a spatula and a pair of tweezers, the brain was gently pushed to ventrally cut the optic nerves to rostrally free the brain and to reveal and cut the pituitary stem. Caudally cranial nerves and meningeal flanges were severed. As a result, the brains shed and fell off under their own weight. The nerves were cut with a vise, and the meninges .
were removed from the pituitary stalk to facilitate infiltration of the fixative solution. The time elapsed from euthanasia to brain and pituitary fixation was approximately 5 min.
. . Histology . . . Fixation and sectioning Each brain was placed in a biosafety container with 1500 ml of 4% paraformaldehyde in 0.1 M phosphate buffer, pH 7.4 at 4 • C. After 48 h, to obtain standardized anatomical planes, the brains were coronally cut from the optic chiasm to 1 cm caudally to the mamillary body, yielding a squared block containing almost the entire diencephalon. Block tissues were fixed for one more week after renewing the fixation solution every 48 h and then immersed in solutions with increasing sucrose concentrations (from 20 to 100%) until they sank. Subsequently, diencephalic blocks were carved to obtain standardized cutting planes and sectioned serially in the coronal, sagittal or horizontal plane at 40 µm using a freezing microtome. For coronal sections, whole diencephalons (left and right sides: four animals) were used, whereas the rest of the diencephalons studied (two animals) were cut through the midline, using each side to prepare sagittal or horizontal sections. In addition, 3600 Nisslstained coronal histological sections were obtained from the four diencephalons, whereas 900 sagittal and 900 horizontal Nisslstained histological sections were obtained from the other two diencephalons. In total, 1080 coronal histological sections were used for calcium-binding proteins (CBPs) immunohistochemistry and 720 coronal histological sections for ionized calcium binding adaptor molecule 1 (IBA1) and glial fibrillary acid protein (GFAP) immunohistochemistry. After fixation, a 3% retraction occurred in the diencephalic blocks. After sectioning, in Nissl-stained sections, the thickness retraction was 4 microns, as calculated by three experienced microscopists using a x100 Leica plan apochromatic objective. Thus, histological sections stained with Nissl had an estimated final thickness of 36 µm.

. . . Nissl staining
Sections were stained by immersion in 1% cresyl violet (C5042, Sigma-Aldrich, Merck KGaA, Darmstadt, Germany) for 10 min at pH 3.0 and in 96% alcohol + acetic acid for staining differentiation. Lastly, the sections were dehydrated in increasing alcohol concentrations, from 50 to 100%, followed by clearing in xylene (3 × 3 min), mounting and cover slipping. The nomenclature used in this paper followed that reported in the Paxinos' atlas of the mouse brain (41). SFR nuclei of this atlas were clearly recognized in our material. In addition, for an appropriate localization of SFR nuclei in the ewe brain, we followed the guidelines based on a recently updated neuromeric model developed in rodents (42) (Figure 1).

. . . Immunostaining
Alternate coronal serial sections were stained for calbindin, calretinin, parvalbumin, IBA1 and GFAP (for details about the antibodies used, see Table 1). Free-floating sections were  (Table 1), for 48 h at 4 • C. Non-specific labeling was blocked using fetal calf serum (10%). After washing three times in TBS-Tx, for 15 min, all sections were incubated with an anti-rabbit biotinylated secondary antibody (biotinylated anti-rabbit IgG H+L, BA-1000; Vector, Burlingame, CA, USA) or with an anti-mouse biotinylated secondary antibody (biotinylated anti rabbit IgG H+L, BA-2000; Vector, Burlingame, CA, USA) at a 1:200 dilution in TBS-Tx for 120 min at room temperature. The sections were then washed with TBS-Tx and incubated for 180 min in avidin/biotinperoxidase (ABC complex, Vectastain Standard ABC kit PK-4000; Vector, Burlingame, CA, USA) and further washed with TBS-Tx, followed by Tris-HCl, pH 8.0. Lastly, the sections were incubated in 3,3-diaminobenzidine tetrahydrochloride (DAB; D-9015; Sigma-Aldrich, St. Louis, MO, USA) with 0.006% H 2 O 2 to visualize the peroxidase reaction. Negative controls, processed without the corresponding primary antibody, were performed to confirm immunostaining specificity; in all cases, the results showed the specificity of the primary antibodies used in this research. Calcium-binding protein cell phenotype analysis across a high number of representative mammalian species showed that the types of positive neurons generally coincided among species. Previous analysis indicated that these proteins have a highly phylogenetically conserved molecular structure in the brain in general and in the cerebellum in particular (30,43,44). Thus, superior colliculus and cerebellum (Supplementary Figure 1) have been analyzed as positive controls, showing immunostaining properties similar to those reported for humans and other high mammals, such as the common marmoset (29). Antibodies used by us are not commercially validated for sheep. However, previous literature . /fvets. .

. . . Image analysis
Panoramic mosaics of entire sections of the rostral half of the brain were captured at selected interaural coordinates. Digital photomicrographs (mosaics) were taken under a Leica DMRB microscope with X10 objectives (Leica Plan Apo), assembled using the "Virtual slice" module of Neurolucida 8.0 © (MBF-Bioscience, Williston, Vermont, USA), and adjusting the microscope illumination source before each image capture. To measure the area and perimeter of the cells, photomicrographs of Nissl-stained sections were processed in ImageJ 2.0 software. Using the automatic method of analysis by density threshold segmentation, a detection range was applied to select the appropriate outlines of cell bodies, which were visually inspected to avoid counting superimposed or confluent groups of cells. To enhance photographs by pseudo-color transformation, improving the delimitation of Nissl-stained panoramic neuronal nuclei, we used different Image J software plugins (ICA 3, Spectrum, and 16-color tables). To numerically assess differences in the density of coronal sections immunostained with IBA1, mean gray levels were analyzed (ImageJ program) in selected frames taken on the horizontal axis at different distances from the 3 rd ventricle.

. . Serial Nissl staining sections of the diencephalic and hypothalamic SFR nuclei
To set cytoarchitectonic boundaries in panoramic highresolution images, nuclei were delimited by exploring differences in cell density, cellular shape and size, and Nissl staining intensity. In this study, the preoptic area (POA), the anterior and the posterior areas of the terminal hypothalamus (THy), were analyzed to more easily describe the localization of SFR (Figure 1). This approach helped us to locate, in a single sagittal section, the four nuclei of interest, ARC in the posterior (Figure 1 dot 1

. . . Intergeniculate leaflet of the thalamus
The thalamic lateral geniculate nucleus (LGN) was large (approximately 1 cm in the main diameter), oval, and located laterally and dorsally in the brain (Figure 2A). Pseudo-color images allowed us to define IGL borders easily and clearly in Nisslstained sections after color conversion ( Figure 2B bottom). Thus, in both sagittal ( Figure 2A) and coronal planes ( Figures 3A, G), IGL was identified by the snake plot shape of a highdensity and hyperchromatic band of cells. This morphological feature facilitated the identification of the thalamic IGL in the border between the dorsal and ventral LGN. However, due to the presence of separate clusters and to its convolved distribution, IGL was difficult to delimit in horizontal sections (not shown). IGL neuronal bodies were identified by their fusiform shape in the coronal plane and by their mediumto-large size (Area M: 649 µm 2 -SD ± 346.6) ( Figure 3G; Supplementary Table 1).

. . . Arcuate nucleus
In sagittal sections, from its rostral limit with the mammillary body, ARC appeared as a ribbon of cells, distributed, along the entire ventricle extension to the infundibulum (Inf) and the median stalk eminence (Figures 2D, G, dotted line). In the ventral THy, in the coronal plane, ARC wrapped around the . /fvets. . of the 3 rd ventricle (Figures 4A-C). In coronal sections, ARC was triangular, extended for approximately 3 mm dorsally, and showed strongly stained neurons ( Figure 3H). Its neurons were medium to large (averaged area: 590 µm 2 -SD ± 290.5), hyperchromatic, triangular and more densely packed in the ventral part (Supplementary Table 1).

. . . Suprachiasmatic nucleus
SCh was located medially and ventrally on the floor of the 3 rd ventricle, bordering the optic tract and chiasma. SCh was small, round, or oval at the horizontal plane ( Figure 4C), and conical at the sagittal ( Figure 2E) and coronal ( Figures 3E-J) planes. According to Nissl staining, SCh can be subdivided in ventral and dorsal areas. SCh neurons had different shapes of neuronal somata, mostly with stellate and oval contours, with small or medium sizes (Averaged area: 161.9 µm 2 -SD ± 181.2) ( Figure 3J;   Figure . (D-F) Three representative sections from the rostral hypothalamus (blue area in Figure ). The labels in the bottom of each image (C to C ) match the level of sectioning corresponding with vertical red lines of Figure .  and oval in coronal sections (Figures 3E, F). Given its large size, this nucleus was observed along more than 20 rostral sections ( Figures 3E, F). Because SO is close to the optic tract, at the horizontal plane, SO was identified at the lateral part of the most ventral horizontal sections ( Figure 4A) and medially to the dorsal ones ( Figure 4C). Thus, the optic tract and optic chiasm border the SO nucleus (Figures 2A, 3E, F, 4A-C), whose neurons are uniformly round or fusiform and large (Averaged area: 915 µm 2 -SD ± 327.6) ( Figure 2K; Supplementary Table 1).

. . Neuronal morphological features and CBPs immunostaining of the diencephalic and hypothalamic SFR nuclei
The neurons of SFR nuclei showed Cb and Cr, but not Pv, immunoreactivity. All CBPs immunocytochemistry results from serial sections in the coronal plane were compared with the corresponding Nissl-stained alternate sections in the coronal plane.

. . . Intergeniculate leaflet of the thalamus
Cb and Cr -stained sections showed, in IGL, restricted immunopositively terminal buttons, which made it possible to delineate the nucleus from LG subdivisions (Figures 5A-D). By immunoreactivity against Cr, dense oval cell bodies were observed (Figures 5C, D), as in Nissl-stained sections ( Figure 3G). In Pv immunostained sections, the intense positive reaction along LGN blurred the nucleus, making it difficult to determine a specific reactivity in cell bodies (not shown).

. . . Arcuate nucleus
In ARC, stellates cells were immunoreactivity against Cb and Cr ( Figures 5E-G). Those showing immunoreactivity to Cb were a specific group of stellate multipolar medium or small neurons with a large number of dendrites ( Figures 5E, F). Cr-immunoreactive neurons were dispersed in the most ventral part, near the third ventricle ( Figure 5G), such as small neurons observed in Nisslstained sections ( Figure 3H).

. . . Paraventricular nucleus
PVN neurons were not immunoreactive to Pv and Cr (data not shown) but were immunoreactive to Cb, which allowed us to observe large multipolar neurons with well-defined cell bodies and well-defined dendrites ( Figures 5H-K). Similarly, in Nissl-stained sections neurons from PVN were large and round ( Figure 3I).

. . . Suprachiasmatic nucleus
SCh contained Cb-immunoreactive neurons and showed positive fibers from the optic chiasm entering the nucleus, generating a dense terminal field (Figures 5L, M). After comparing immunoreactive against Cb with equivalent Nissl-stained sections, only a few of the Nissl-stained neurons ( Figure 3J) were positive for this antibody ( Figure 5L).

. . . Supraoptic nucleus
In our material, most SO neurons were immunopositive again Cb and/or Cr (Figures 5N-P). Immunopositively stained dendrites and cell bodies showed strong Cb and Cr staining (Figures 5N-P). Its medium-to-large neurons shows strongly immunopositive dendrites. Immunostained sections showed much smaller number of cells than the Nissl-stained ones (compare Figures 5N-P with Figure 2K).

. . Glial architecture IBA and GFAP immunohistochemistry
No distinct features of astroglia and microglial cells were detected in IGL with respect to the remaining lateral geniculate body (not shown). In panoramic views, glial cell size and cell density were increased in extensive areas around ventricles in POA (rostral to the hypothalamus) and PHy, both for IBA 1 and .
/fvets. . GFAP immunocytochemistry (Figures 6, 7). In our material, medial (MPA) and lateral (LPO) regions of POA were distinguished based on differences in staining. IBA 1 immunoreaction labeled both areas, but GFAP immunoreaction only labeled MPA (compare Figures 6, 7). Inside the POA and PHy, an intense immunoreaction in microglial and astroglial cells was observed inside SFR nuclei ( Figures 5, 6 dotted lines). In PHy, a high concentration of IBA1-immunoreactive cells was identified along the border of the 3 rd ventricle to the ventral area around and inside ARC ( Figures 6A, B dotted line); however, the density of immunoreactive cells decreased in the infundibular area, allowing us to define an ARC microglial-specific region separate from the periventricular zone ( Figures 6A, B).
In IBA1 immunostained sections, from PVN and the 3 rd ventricle, to the lateral limit of LPO, cell density decreased gradually dorsally and laterally, as shown after comparing gray density values from selected square areas in Figure 5C  and dorsally slight decrease in immunoreactive microglial cells ( Figure 6).
GFAP immunoreactivity also showed a high concentration of astrocytes around the ventricle and POA, but the densest reactive areas were more restricted to the paraventricular zone ( Figure 6C). As in IBA1 sections, GFAP immunostaining density decreased laterally in POA, matching MPA and LPA borders (e.g., please note these differences in cell density in details 2, 3 and 4 in Figure 7). In GFAP-immunostained sections, PVN neurons appeared as negative dots framed by densely stained astrocytic feeds ( Figure 7C square 3). In the most caudal, sections immunostained against GFP show tanycytes which cross the stalk and the median eminence forming a superficial well-stained limit and a superficial barrier ( Figures 7A,  B squares 1 and 2).

. . Macroscopic localization of nuclei in Nissl-stained sections
Measurements of the main central axis of the mammillary body in the three sectioning planes in histological Nissl-stained sections, compared with measures in macroscopic images, show small differences in retraction (<3%). As mentioned above, with measurements of the thickness of 20 randomly selected sections taken by three different observers (z-axis, and using a Leica X100 plan apo objective), we calculated a shrinkage of 4 (+/-) µm (final averaged thickness 36 µm). These data facilitated slight corrections in the translation of coordinates using the fixation and sectioning protocol applied here. Moreover, when coordinates are translated from frozen specimens, which is the most suitable approach to tissue extraction for molecular analysis, the measurements must .
/fvets. . be modified after evaluating changes in tissue contraction in each procedure. For micro-dissection, coordinates for extraction were defined in present paper in a caudal-to-rostral direction taking as referent point the anterior limit of the mammillary body (easily located macroscopically by its spherical shape close to the midline). From this reference point, ARC was localized in a sagittal view at 1.5 mm at the vertical axis, 2 mm at the horizontal axis ( Figure 8A), and 1 mm laterally to midline ( Figure 8B). By using as a reference, a parallel line to the ventral surface (tangential to the ventral most point of the mammillary body and the surface of the optic tracts -green line in Figure 8A), the periventricular region was located at its main axis oriented in an angle of 47.43 o . By using the same reference line, SCh specimens were localized 0.7 mm dorsally to the recess ( Figure 8A, red arrows), and PVN was dissected by cutting a 4.2 × 1 mm rectangle ( Figure 8A rectangle). Due to its localization, sandwiched between the optic tract (OT) and the preoptic area, SO coordinates cannot be accurately defined, but a tentative localization is suggested in Figure 8B (black dotted line). Furthermore, by its convoluted structure and oblique orientation, IGL coordinates are very difficult to precisely define in macroscopic views; nevertheless, a tentative localization is also provided in Figure 8 (red dotted rectangle).

. . General considerations
After SFR nuclei activation by changes in photoperiod, peptide regulation induces neuronal activation and plasticity (11)(12)(13) leading to morphological changes in the hypothalamus. These changes in the estrous cycle, previously studied in rats (49) and as a function of seasonal fluctuations, mainly regulated by melatonin (33) may induce size and shape differences in SFR nuclei. In this study, ewes were euthanized under specific environmental, nutritional, and estrous cycle conditions. However, potential changes in morphology may be induced by plastic reorganization after changes in hormonal induction. Therefore, the anatomical data reported here must be revised under other estral conditions. Relative SCh, ARC and PVN differences in shape and localization were identified in this study. SO showed substantial qualitative differences in both its layout and relative size in ewes, compared with other mammals (e.g., rodents, primates. pigs) (26, 31); however, the basic neuroanatomical plan of organization of ewe SFR nuclei, is similar to that of other mammals.
In line with differences in reproductive functions between rams (e.g., mounting) and Ewes (e.g., pregnancy and lactation), strong anatomical differences have been previously shown in the nuclei studied here (sexual dimorphism) (50, 51). Accordingly, in the future a specific anatomical systematization should be performed in rams for a complete neuroanatomical description of nuclei involved in seasonality in sheep (Ovies orientalis aries).
KNDy neurons form an interconnected network that fires synchronously to drive GnRH release during a pulse, inducing combined differences between nuclei closely positioned in the hypothalamus (23,24). For future molecular analysis, functional and anatomical interconnections between hypothalamic nuclei involved in seasonal sexual activation should be studied using welldefined samples of tissue restricted to specific nuclei. A study using nuclei microdissection to analyze large areas of the hypothalamus quickly immersed in liquid nitrogen reported excellent results about changes in GnRH-R mRNA levels after stress in merino sheep (52). Thus, in the present study, we provide a cytoarchitectonic guide for a more precise microdissection of SFR nuclei, enabling researchers to collect more appropriate samples for conventional proteomic or genomic techniques, individual cell sequencing or chemical analysis of hormones, transcription factors, and peptides. Preparing serial sagittal sections after micro-dissection enables to confirm the appropriate extracted areas, as previously we reported in the rat (53).

. . Intergeniculate leaflet nucleus
IGL was first described in rats (54), and golden hamsters (55) as a well-defined band of cells located between the dorsal and ventral medial geniculate bodies. When comparing rodents and higher mammals, such as primates (26) or sheep (present results), the well-defined straight band between dorsal and ventral subdivisions of the visual thalamus were masked and hidden by the complex laminar organization of the dorsal geniculate bodies. However, in our Nissl-stained sheep sections (Figure 2A), large fusiform neurons, delimiting IGL, were more clearly observed in the coronal plane ( Figure 3G). In rodents, nevertheless, its architecture is well defined, but both the localization and limits are difficult to distinguish in some gyrencephalic vertebrate species (26) because primates once had a generic wide pregeniculate nucleus (PGN) included both the ventral geniculate nucleus and IGL (56). In the visual thalamus of the rock cavy (Kerodon rupestris; rodent), IGL is identified as a well-defined and straight band between dorsal and ventral subdivisions, but in marmosets (primates), this subdivision is masked and hidden by the complex laminar organization of the dorsal part of the LG (57). Due to its intricate architectonic organization, the primate IGL should be neuroanatomically better defined by using additional neurochemical markers (e.g., VIP immunocytochemistry). However, in the ewe, positive Cb and Cr terminal fields, most likely coming from the retina" allowed us to trace the borders of this subdivision. In summary, due to its larger fusiform neurons and Cb and Cr immunoreactivity ( Figures 5B, C), IGL can be well differentiated from the remaining visual thalamus in the ewe ( Figure 3A).

. . Arcuate nucleus
Anatomically, six types of neurons have been described in previous studies using the Golgi method in the rat ARC nucleus (58). Different neurochemical phenotypes have also been reported, particularly a neuronal type in which kisspeptin, neurokinin A and dynorphin coexist, known as KNDy neurons (59). These key neurons regulate seasonality outputs for reproduction (kisspeptin, neurokinin, dynorphin) by activating GnRH pulses in sheep (32) Because these neurons are mainly located in ARC (23, 60), techniques demonstrating kisspeptin, neurokinin or dynorphin (not available in our material) should be used as complementary markers to be able to better correlate the morphology and topography of this nuclei with present results. In ARC, several chemical neuronal types have also been described; thus, an accurate analysis in the ewe will require a peptide immunocytochemistry analysis, which was not performed in our study. Nevertheless, in our material, the morphology (Nissl) and immunoreactivity (CBPs) of dorsal ARC neurons showed that they were roughly similar in size and location to KNDy neurons.
. . Paraventricular nucleus PVN shows the highest peptidergic and glial activity in the hypothalamus and is located at the border of the third ventricle.
Surrounding this structure, some regions of special interest for their association with seasonal reproduction in sheep, such as the set of A15 dopaminergic cells, regulate the negative feedback of estradiol, which in turn is another compound that regulates both GnRH production and pulses (32, 61, 62). In ewes, the columnar distribution at the sagittal plane allowed us to distinguish strongly Nissl-stained neurons into three potential subdivisions ( Figure 2E).

. . Suprachiasmatic nucleus
Light cycle-sensitive pacemaker SCh neurons are directly regulated by retinal and indirectly by IGL projections through the geniculate thalamic tract (63). SCh, a highly phylogenetically conserved structure with a similar location in many species of vertebrates, is placed between the 3 rd ventricle and the optic chiasm and contains small and packed neurons (64). In ewes, we also observed small neurons (∼430 µm 2 ) and strong immunoreactivity for fibers containing Cr or Cb presumably belonging to retinal projections ( Figures 5L, M).

. . Supraoptic nucleus
SO interconnects several of the seasonality nuclei (PVN and ARC), directly projecting to the pituitary gland and to the region surrounding the third ventricle (44,45) Because SO is involved in glial activation, we will analyze its morphological features more extensively in the next section (51). Despite differences in animal size, SO nucleus and its neurons are much larger in ewes than in rodents, but they have the same topographical distribution (26, 41).

. . Glial architecture
To our best knowledge, this is the first analysis of glial cell architecture in the ewe hypothalamus and diencephalon. Glial cells are widely known to play a key role in the synthesis, recycling and delivery of hypothalamic hormones and factors (65) and thus GnRH in the hypothalamus of different mammals, including sheep (51,66,67). Increases in cell density and cell size around the 3 rd ventricle shown in our material after IBA1 and GFAP immunoreactivity, indeed expresses glial overactivation (increase cell size, number of profiles, immunostaining and cell density) which can be related to potential areas of higher metabolic regulation of sexual hormones (GnRH) (68). Accordingly, the densest areas of microglial and astroglia immunoreaction around PVN and SO, in the anterior hypothalamus, and around the ARC, in the THy, may be related with a heightened endocrine or paracrine regulation.
Notwithstanding the above, both types of reactive glial cells (as shown by larger and denser population of microgliocytes and astrocytes) are found around neurons. These findings suggest the potential role of denser immunoreactive areas with overactivated paracrine secretion. When comparing IBA1 and .
/fvets. . GFAP immunoreactivity, microglial reactivity was more extended in all section levels, suggesting a microglial more widespread in areas potentially implicated in secretory regulation (please compare Figures 6, 7). Tanycytes formed a dense ventricular barrier of the 3 rd ventricle in the stalk-median eminence ( Figures 7A, B squares 1, 2). Due to the absence of any other reactive glial cells, this finding suggests that these areas, as in other species, regulate the access of metabolic signals to the hypothalamus (69).

. . Concluding remarks
Our results establish a guide for future histological, physiological and molecular experiments related to hypothalamic effects of interventions in farm animals aimed at improving fertility regulation (e.g., GnRH neurons, KNDy cells, and melatonin receptors). The data provided here may also have translational relevance because the sheep brain, considering its similarity in anatomy and organization, is a good experimental model for human neuroanatomical research, as recently suggested (70). Finally, by establishing a procedure for the section-to-brain translation of coordinates, this research may also help us to develop suitable procedures for nuclei/tissue extraction from specific nuclei of interest toward assessing changes in mRNA levels by RT-PCR or in protein levels by Western blot.

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement
The animal study was reviewed and approved by the University of Salamanca.

Author contributions
MM and CP designed the experiments. MM and IP perform the experiments. MM and RC wrote the paper. CP and JA participated in the discussion of the results and corrected the manuscript. All authors contributed to the article and approved the submitted version.

Funding
This study was supported by the Spanish Ministry of Science and Innovation (Ministerio de Ciencia e Innovación -MICINN) Grant number: PID2020-117266RB-C21.